Self-organized criticality in an interface-growth model with quenched randomness 
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We study a modified model of the Kardar-Paris- Zhang equation with quenched disorder, in which 
f"^ , the driving force decreases as the interface rises up. A critical state is self-organized, and the 

' anomalous scaling law with roughness exponent a ~ 0.63 is numerically obtained. 
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Nonequilibrium interface growth and its scaling properties have intensively been studied. [l[ There are a number 
of applications of random interface-growth problems in molecular beam epitaxy, bacteria-colony growth, and fluid 
' invasion in porous media. In many rough interfaces, the root-mean-square (rms) width W(l,t) obeys a dynamic 
\ scaling law: 

,— , w(i, t) = ((h( x , t) - (h( X , t))) 2 ) 1 ' 2 ~ i a h(t/i% (i) 

O ' where h(x,t) is the surface height at time t, z = a//3 is a dynamic exponent, / is an interval size, and the scaling 

S function h{u) satisfies asymptotically h(u) ~ for u <C 1 and h{u) -> 1 for ti > 1. The exponent a is called the 
roughness exponent and (3 is the growth exponent. It has been widely believed that many systems lie in the same 
+L> . universality class as the Kardar-Parisi-Zhang (KPZ) equation. [2[ The (1 + 1) KPZ equation has a form: 
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where v and A are constants, and T)(x,t) denotes the Gaussian white noise satisfying (77} = and (rj(x,t)ri(x' ,t')) = 
1 ■ 2DS(x — x')5(t — t'). The roughness exponent a is 1/2 and the growth exponent 8 = 1/3 in the (1 + 1) KPZ 
equation. However, many experiments give larger values of a in the range 0.6-0.9 @-[6i- To explain the anomalous 
exponent, some authors studied modified models of the KPZ equation, in which the noise was assumed to have a 
power-law distribution, 0] although the physical origin of such a noise remains unclear. The randomness is quenched 
in space in many experiments, Thus, the KPZ equation with quenched disorder (KPZQ equation) was also studied 
theoretically. The KPZQ equation in one dimension is written as 
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^j! . where r)(x,h) denotes the quenched disorder satisfying (rj) — and (r/(x, h)rj(x' , h')) — 2DS(x — x')S(h — h'). The 
nonlinear coefficient A is proportional to the average velocity v from the kinematic origin in most case, however, such 
a nonlinear term appears in anisotropic media even for v = 0. jlOl ] 

The KPZQ equation exhibits a pinning transition when the driving force / is decreased. When the driving force 
/ is larger than a critical value f c , the interface grows with an average velocity. When / is below / c , the interface 
is pinned in the quenched random medium. At the critical driving force f c , the scaling law W(l) ~ l a is satisfied, 
and the exponent a is evaluated at a ~ 0.63. The exponent a is related to two exponents and v± in the directed 
percolation problem as a = v\\/v_i_, where v\\ is the exponent for the longitudinal correlation length £1 and u± is the 
exponent for the transverse correlation length (ll| The exponent a — 0.63 is close to the experimental exponents 
in paper wetting [ilj . The scaling law is theoretically satisfied only at the critical driving force f c . However, the 
scaling law is observed even if the parameter is not precisely controlled to be the critical value. The critical condition 
might be to be self-organized by a certain unclear mechanism. An interface-growth model which exhibits the self- 
organized criticality (SOC) was proposed by Sneppen [12} . In his algorithm (the model B in [12]), the site with the 
smallest random force r/(x, h) is selected , and the updating h —± h + 1 is done at the site. Then, neighboring sites 
are sequentially updated as h — > h + 1, if \h(x) — h(x — 1)| < 1 is satisfied. Such a chain reaction rule is used in many 
models exhibiting the self-organized criticality. [l3| A self-organized random interface with a — 0.63 was obtained in 
the model of Sneppen, but the physical origin of the updating rule is not clear. Thus, we think that the problem of 
the anomalous exponent in random interfaces is not completely solved yet. In this brief report, we propose another 
simple mechanism of the self-organized criticality in a modified model of the KPZQ equation. 

We consider an experiment of paper wetting as performed in [11]. The driving force in this problem is the surface 
tension F . The paper is dipped into a basin filled with suspensions of ink or coffee. When the interface rises up, 
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FIG. 1: Critical driving force f c (solid line with rhombi) in the discrete KPZQ equation (6) at v = 5 and L = 10000 for 
A = 0, 0.25, 0.5, 0.75, 1, 1.25 and 1.5. The marks + denote the average values of F — d(h) for the same values of A in the time 
evolution of Eq. (7). 

the upward driving force decreases owing to the gravity as / = F — d$h, where do is proportional to the density p 
multiplied by the gravitational acceleration g, because the gravitational potential increases with the square of the 
interface height. [141 ] . Furthermore, the evaporation from the surface of the paper needs to be taken into consideration, 
if the humidity is low. As a rough approximation, the interface height decreases on the average by —eh, because the 
volume of the liquid is proportional to hLw where w is the average thickness and L is the side length of the paper, the 
evaporation rate is proportional to the area 2hL of two surfaces of the paper, and the evaporation rate is proportional 
to the decrease rate of the liquid volume. That is, the evaporation effect also decreases the upward driving force as 
/ = F — eh as a first order approximation. Thus, our model equation is written as 



where d = do + e. If the term dh is approximated at the spatial average d{h), where (h(t)) is the spatial average of 
the interface height in system size L: (h(t)) = (l/L) j Q L h(x,t)dx, our model equation is expressed as 



This model equation is closely related to the KPZQ equation, or it can be interpreted as a modified KPZQ equation 
with a feedback term. If F > f c , the interface grows upward from the flat initial condition h(x,0) = 0. The driving 
force becomes weaker, as the average height (h) increases. The interface growth is expected to stop when F — d(h) ~ f c 
is satisfied owing to the pinning transition in the KPZQ equation. After the interface is pinned, the critical condition 
F — d(h) ~ f c is maintained forever, because (h) is constant in time. Thus, the self-organized criticality can be realized 
by this very simple mechanism. 

We have performed direct numerical simulations of Eq. (3) and (5) using a simple Euler method. The KPZQ 
equation (3) is discretized as 



h(x, t+At) = h(x, t)+At [v{h(x + 1, t) - 2h(x, t) + h(x - 1, t)} + \{h(x + 1, t) - h(x - 1, <)} 2 /4 + / + rj(x, h)] . (6) 
The modified model Eq. (5) is discretized as 

h(x, t+At) = h(x, t)+At [v{h{x + l,t) - 2h(x, t) + h(x - 1, t)} + \{h{x + l,t)-h(x- 1, i)} 2 /4 + F — d(h) + r/(x, h)] . 



The parameters are fixed to be v = 5 and At — 0.01. We have used quenched disorder r)(x, h) which is uniformly 
distributed between and 3. @ Figure 1 shows a critical curve f c (solid line) of the discrete KPZQ equation (6) 
as a function of A. The system size is L = 10000 and periodic boundary conditions h(L,t) = h(0,t) are assumed in 
the numerical simulation. The pinning phenomenon is observed below the critical curve. The critical driving force f c 
decreases with A. 

Figure 2(a) shows the time evolution of h(x, t) for L = 10000, A = 1, d = 0.00005 and F — 0.17 by Eq. (7). Snapshot 
profiles of h(x, t) at 1000 X n (n = 1, 2, • • • , 12) are shown in Fig. 2(a). The initial condition is h(x, 0) = 0. For small 
t, the interface rises up randomly. As the interface rises up, the driving force becomes weak and finally the interface 
is pinned. Near the pinning transition, the interface growth occurs intermittently both in space and time, that is, the 
interface is pinned in most positions and the interface growth occurs locally. Figure. 2(b) shows the time evolution 
of F — d(h) for the numerical simulation shown in Fig. 2(a). The driving force F — d(h) decreases and finally takes a 
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FIG. 2: (a) Successive snapshot profiles of h{x, t) at 1000 x n (n = 1, 2, ■ • ■ , 12) by Eq. (7) at v = 5, A = 1, F = 0.17, d = 0.00005 
and L = 10000. (b) Time evolution of F — d(h) for the numerical simulation shown in (a). 
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FIG. 3: Relation of W(L) as a function of L after the interfaces are pinned for Eq. (7) at v = 5, A = 1, F = 0.17 and d — 0.00005. 
The dashed line is the fitting curve W(L) = 0.0907L ' 642 by the least squares method. 

constant value by the pinning. The final value is 0.157 in this quenched disorder. The final values slightly depend on 
the quenched disorder. We have calculated the average value using 50 random quenched media. The average values 
of F — d(h) at the stationary state for various values of A are plotted by + marks in Fig. 1. The stationary values 
of F — d(h) almost locate on the critical curve / C (A) of the KPZQ equation. This implies that the self-organized 
criticality is realized in our model. We have calculated numerically the rms width W(L) after the interfaces are 
pinned for various values of L at A = 1, F = 0.17 and d = 0.00005. Figure 3 shows a relation of W(L) as a function 
of the system size L. The exponent a is evaluated by the least squares method as a = 0.64, which is close to 0.63 for 
the KPZQ equation at the critical point. The anomalous critical exponent is realized at the stationary state in our 
model, even if the initial driving force F is not precisely controlled to be f c . 

To check the validity of the mean- field approximation (5) to Eq. (4), we have performed direct numerical simulation 
of the discretized model of Eq. (4) with d = 0.00005. Figure 4(a) shows the time evolution of F — d(h) at v — 5, A = 
l,F = 0.17 and L = 10000. The time evolution is very close to the one shown in Fig. 2(b). Figure 4(b) shows the 
relation of W(L) vs. L. The exponent a is evaluated by the least squares method as a = 0.64. These numerical 
results show that the approximation by Eq. (5) is rather good in our simulation at sufficiently small d. However, 
it is expected that the scaling properties on very large length scales in Eq. (4) are different from the ones with the 
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FIG. 4: (a) Time evolution of F-d(h) by the discretized model of Eq. (4) at v = 5, A = 1, F = 0.17, d = 0.00005 and L = 10000. 
(b) Relation of W(L) as a function of L after the interfaces are pinned for Eq. (4) at v = 5, A = 1, F — 0.17 and d = 0.00005. 
The dashed line is the fitting curve W(L) = 0.0917L 0,639 by the least squares method 
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roughness exponent a = 0.63. The crossover length is expected to be 0( y/v/d), above which the linear damping term 
— dh in Eq. (4) becomes dominant compared to the diffusion term vd 2 h/ dx 2 . On the other hand, in the mean-field 
model (5), the damping term acts only for the spatial average (h), and does not have a direct effect for the spatial 
fluctuation of h(x,t). 

To summarize, we have found a very simple mechanism of self-organized criticality in a modified model of the KPZQ 
equation. The critical pinning state is naturally realized by decreasing the driving force. Both the gravitational effect 
and the evaporation effect arc considered to decrease of the driving force effectively The gravitational effect is 
controlled by changing the inclination angle of the paper and the evaporation effect is controlled by the humidity. 
This might be a mechanism of the anomalous scaling in paper wetting experiments. Similar mechanisms might work 
even for other experiments of interface growth in quenched random media. For example, there is another feedback 
mechanism for the self-organized criticality in the paper wetting problem where the paper is set horizontally and 
the gravitational effect does not work, and the evaporation effect is neglected. The average growth velocity v of the 
interface decreases as £ -1 / 2 (Washburn's law) because of the viscous resistance force in proportion to (h). [H|, EH 
Then, the A parameter decreases effectively in time owing to the kinematic effect, and the self-organized critical pinned 
state will be realized, because the critical driving force f c is a decreasing function of A in the KPZQ equation as shown 
in Fig. f . Furthermore, in general, there might exist other systems, in which critical states are self-organized by some 
feedback terms. We would like to study such self-organized criticality by the feedback effect in the future. 
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